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Abstract 

The two-body Dirac equations for the bound qq systems are obtained from the dif- 
ferent (five) versions of the 3D-equations derived from Bethe-Salpeter equation with the 
instantaneous kernel in the momentum space using the additional approximations. There 
are formulated the normalization conditions for the wave functions satisfying the obtained 
two-body Dirac equations. The spin structure of the confining qq interaction potential is 
taken in the form xj® <8> 72 + (1 — x)I\ (8> 1*2, with < x < 1. It is shown that the two-body 
Dirac equations obtained from the Salpeter equation does not depend on x. As to other 
four versions such a dependence is left. For the systems (its), (cit), (cs) the dependence of 
the stable solutions of the Dirac equations obtained in the different version on the mixture 
parameter x is investigated and results is compared with such dependence of 3D-equations 
derived from Bethe-Salpeter equations without the additional approximation and some 
new conclusions are obtained. 
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1 Introduction 



The Bethe-Salpeter (BS) equation provides a natural basis for the relativistic treatment of 
bound qq systems in the framework of the constituent quark model. But due to fact that 
the BS wave function (amplitude) has not probability interpretation, three-dimensional (3D) 
reduction is necessary. Review of investigations of bound qq systems (mesons) on basis of 
equations obtained in different versions of 3D-reduction of BS equation in the instantaneous 
(static) approximation for kernel of BS equation is given in Ref. [1]. In literature there are 
known five such versions formulated in Refs. [l]-[7], below noted as SAL [2], GR [3], MW [4], 
CJ [5] and MNK[6], [7] versions. The last four 3D-equations have correct one-body limit (the 
Dirac equation) when the mass of one of the particles tends to infinity. As it is well-known the 
Salpeter has not such a limit. Note that Gross equation is obtained only for m>\ ^ m 2 case, 
while other versions work for the equal masses (mi = m 2 ) too. 

Below we shall consider the problem how to get two-body Dirac equations for the bound qq 
systems from above mentioned 3D-relativistic equations, and how to formulate corresponding 
normalization conditions for the wave functions. Then, these equations will be used for the 
investigation of some aspects of the problem connected to the mass spectra for bound qq 
systems (mesons). Namely, the dependence of the existence of stable solutions of obtained 
equations and of the mass spectra on the Lorentz(spin) structure of confinement qq interaction 
potential will be studied. Further, the comparison to the results obtained without additional 
approximations will be made. 



2 The two-body Dirac equation for bound qq systems 
and normalization conditions for the corresponding 
wave functions 

To derive such an equation note that all 3D-equations given in Ref. [1] can be written in the 
common form (c.m.f.) 
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Note that from the eq.(l) with the operators II (2) immediately follow the system of equa- 
tions (3.61) in Ref.[l] with definition (3.61-63), if eq.(l) is multiplied from left by projection 
operator A^ 1 " 2 '' and are used their properties: 



A (aia 2 ) = A (ai) ^ ^i) = <^±^ kf i] = S^A™ 

n SAL = a{ + 2 +) - a[v ] = l{t 1 +^)^ GR = 4r } + A<r> = \ (l + £) 

Now if in the operator H SAL use the approximation 
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then we obtain the two-body Dirac equation 
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which already was used for bound qq systems in Refs.[8],[9]. 
In approximation (4) from (2) follows 
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As to MW, CJ and MNK versions for derivation of corresponding two-body Dirac equations 
the additional to (4) approximation is need, namely 



n(M ; 0) n(mi + m 2 ; 0) = n 



(8) 



which is quite natural because it corresponds to zero approximation in iteration procedure for 
solving nonlinear over M eq.(l) for the MW, CJ and MNK versions. As a result from (2) can 
be obtained 
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Thus, we have the following two-body Dirac equations obtained from (1), (2) 
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[M - - h 2 (-p)}V M (p) = n o7 i° ® 7 2 ° 



(ii) 



where the operator n is given by the formulae (6,7,9,10). 

Note that there is the another approach for formulation the two-body Dirac equations, 
namely, generation of the one-body Dirac equation to two-body one, using constrain dynamics 
and relation to quantum filed theory. Review of such an approach is given in Ref.[10]. 

Representing the wave function \I>m(p) as sum of "frequency" components 



*m( P ) = e Air 2) (p)*M( P ) = e ^ 1Q2) (p) 
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from the eq.(ll) follows the system of the equation for the functions 
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Taking the qq interaction operator V in the form [1] (combination of one-gluon exchange 
and confining part of potential) 
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and representing the function \l/^ ia2 ' ) as 
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then for the wave functions Xm" 2 ' from (13) can be obtained the following system of equations 



(16) 
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It is very important that the two-body Dirac equation (17) with effective potential (18) 
with (19) obtained from the equation (1), corresponding to SAL version (2) does not depend on 
parameter x interned in the interaction operator (14), which means that from this equation can 
not be obtained any information on the Lorentz (spin) structure of the confining qq interaction 
potential (14). Second interesting result is that the wave functions satisfying the two-body 
Dirac equation (17) with effective potentials (18) with expression (19, 20) obtained for SAL 
and GR versions (2) of the 3D-relativistic equations have all nonzero "frequency components" 
whereas two components of the wave functions satisfying the equation (1) with projection 
operators (2), are zero, namely: 
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which directly follows (and is well known) from the eq.(l) if it is multiplied (from left) by the 
operators A^ T \ a[ 2 ^ and used the formulae (3). 

For formulation normalization condition for the wave function (12) which satisfies the equa- 
tion (11), we note that normalization condition for Salpeter wave function obtained in Ref.[l] 
(see relation (3.14)) can be written in the form 



< $ m \I1 sal \$m >= 2M 



(25) 



The analogous condition can be derived for wave function satisfying Gross equation (1, 2) 
if we use equation for full Green operator corresponding to the equation (1) 

G = £onr + g UG, g — [M — hi — h,]' 1 , T = 7l ° ® 7°, U = UT V (26) 

Assuming that the operator G~ x exists (being natural at any rate in the bound states, we need) 
from eq.(26) after some transformations can be obtained the following relation 
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Noting that U GR Il GR = Uqr from (27) we have 
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Now using the spectral representation of Green operator G 
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from (28) can be obtained the relation 

< $ M |n G/? |i» M > n Gi? = 2MU GR (30) 

It means that the normalization condition analogous to (25) 

< $A f |n Gi? |$M >= 2M (31) 

holds only in corresponding subspace of the Gilbert space. Note that the condition (31) can 
be obtained from the formula (3.28) of ref.[l], which was not derived, but supposed with an 
analogy to (3.14). 

Now, noting that the two-body Dirac equations (11) for the SAL and GR versions of the 
3D-relativistic equation (1) were obtained in the approximation (4) for the projection operators 
Tl SAL and U GR , the corresponding condition for wave function can be obtained from (25, 31) 
by replacement II s AL =>• Il^ AL (6) and Il GR =>• 11^ (7). Thus we have 



< ^ m \Ho AL ' GR \^m >= 2M (32) 

where Uq AL and U GR are given by formulae (6), (7). 

Further, noting that the projection operator U^f w (9) satisfies condition n^ fH/ n^ fl/,/ = 1, 
from relation (27) can be obtained the normalization condition analogous to (32) i.e. 



< y M \Il™ w \y M >= 2M. (33) 

As to normalization conditions for wave functions satisfying the two-body Dirac equation 
(11), corresponding to the CJ and MNK versions, they can not be derived analogously because 
the corresponding projection operators LTo (10)does not satisfy the conditions Hollo = IIo or 
n n = 1. But bellow we assume (suppose) that the condition analogous to (33) can be written 
in common form 



<^ M \U \^/ M >=2M. (34) 

where operator n is given by the formulae (6, 7, 9, 10) for all versions. As a result with an 
account of the formulae (12, 15, 16) the normalization condition for the components of the 
wave functions Xm 1 "^ takes the form 
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where R^lsj 2 \p) are corresponding radial wave functions. And the potential functions Vog(p, p')> 
Vc(p, p') are represented in form (local potentials) 



(37) 



Now we use the partial-wave expansion for the function Xm i02 \p) [1] 
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j£(x) being the spherical Bessel function. Then from the system of equations (17), the effec- 
tive potentials of which is defined by the formulae (18-23) we obtain the following system of 
equations for the radial functions R^lsT\p) 
SAL version 
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GR version 
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It is interesting to compare the system of equations (41-46) with the system of equations 
obtained from (1) without the approximation (4, 8) (see eqs. (4.16, 17) in [1], neglecting the 
terms corresponding to t'Hooft interaction ) 
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Note that the last expression in (50) is obtained from (3.62) in [1] after some transformation. 

Main difference between the system of equations (41-46) and (48, 49)with the expression 
(50) is following: 1) In the wave functions Rl'sj 2 ^ satisfying the system of equations (48, 49) for 
(SAL) and (GR) versions the nonzero functions are only R^j and respectively (about 

this fact was mentioned above), whereas in the corresponding system of equations (41), (42), 
(43), (44) all components of wave functions Rl'sj 2 ^ are nonzero; 2) The system of equations (48, 
49) for the MW, CJ and MNK versions are nonlinear over M, whereas the system of equations 
(45, 46) are linear one; 3) Dirac equations (41,42) obtained from the Salpeter, equation do not 
depend on x. 
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3 Procedure for solving the obtained equations 



For solving bound-state equations (41-46) or (48,49), we need to specify the interquark inter- 
action potentials Vq-g an d Vc (14). Below for Vc(r) we use the following form [1], [11] 

V c (r) = \a s {m\ 2 ) (—=^^== - V ) (51) 



a s (Q 2 ) 



127T 



33 - 2n f L 



m 1 m 2 , . 

, m 12 = mi + m 2 , H\ 2 = (52) 

m 12 



where Q 2 is the momentum transferred and the |comes from the color- dependent part of the 
qq interaction, n/ is the number of flavors (rif — 3 for u, d, s quarks; rif — 4 for u, d, s, c; rif = 5 
for u, d, s, c, b). uq, Aq, V$ and A are considered to be the free parameters of the model. The 
potential given by expression (51) effectively reduces to the harmonic oscillator potential for 
the light quarks u, d, s and to the linear potential to the heavy c, b quarks if the dimensionless 
parameter A is chosen small enough. Moreover, asymptotically, for a large r it is linear and 
almost flavor-independent. The one-gluon exchange potential is given by standard expression 

[i], ("I 

Voc (r) = -i^SM (53) 

Now we have to specify the numerical procedure for solution of the systems of radial equa- 
tions (41-45), (48), (49). A possible algorithm looks as follows: we choose the known basis 
functions denoted by R n ^{j)). The unknown radial wave functions are expanded in the linear 
combination of the basis functions 



4r 2) (P) = ^2M(2nf £ Ci a L ^ R nL ( P ) (54) 

n=0 

where C^g]^ are the coefficients of the expansion. The integral equation for the radial wave 
functions is then transformed into the system of linear equations for these coefficients. If the 
transaction is carried out the finite system of equations is obtained that can be solved by using 
conventional numerical methods. The convergence of the whole procedure, with more terms 
taken into account in the expansion (54) depend on the successful choice of the basis. In case of 
the confining potential of form (51) it is natural to take as a basis the functions corresponding 
to oscillator potential, which is obtained from (51) at Aq = 0, in non-relativistic limit of the 
system of equations obtained from from (41-46), (48), (49). The radial wave functions in this 
case have the form [l](the formula (4.52)). 



Rul(p) = p 3/2 RnL(z), po = Jni2U>o^a s (ml 2 ), z = £ 



Rul(z) = c nL z L exp{-^) 1 F 1 {-n, L + I z 2 ), c nL = 
where ±Fi denotes the confluent hypergeo metric function. 



(55) 
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Now, satisfying the expression (54) into the system of equations (41-46), (48), (49), the fol- 
lowing algebraic equations for the coefficients C^sj^ can be obtained 

MC^s? = E E flSSWffl (56) 

aia' 2 n'L'S' 

It is necessary to note that the matrix H a p(M) depends on meson mass M only for MW, CJ 
and MNK versions as it can be seen from equations (56) for M is not linear one and therefore 
should be solved, e.g. by iteration. As to the system of Dirac equations (41-46) such a problem 
does not exist. 

4 The numerical results and discussions 

The main problem we have investigated at first stage is dependence of the existence of stable 
solutions of the eq. (56) i.e. the equations (41-46), (48), (49) on Lorentz (spin) structure of 
the confining qq interaction potential, i.e. on the parameter x. This will be done taking as 
examples the us, cu and cs, bound states with constituent quark masses m u = rrid = 280MeV, 
m s = 400MeV, m c = 1470MeV and the free parameters of the confining potential (51,52) - 
luo = 710MeV, V = 525MeV, A = 0.0270, A = 120MeV. 

Note, that in [11] only the SAL version of 3D-reduction of Bethe-Salpeter equation was 
considered as to MW, CJ and MNK without additional approximation (4) with oscillator like 
potential (A = in (51)) were considered in Refs. [7], [12]. 

The results of the calculations are given for states +1 Lj (note, that for cases 3 Si, 3 P2, 3 -Pi 
are neglected additional corresponding terms 3 D±, 3 F 2 , 1 Pi, because they give small contribution 
in the calculated mass). 

The additional conclusions to pure theoretical results formulated at the end of section 2, 
which follow from the tables, are the following: 

• The area of dependence on parameter x existence of stable solutions of corresponding 
equations is a little extended for corresponding Dirac equation. 

• The results for CJ and MNK versions for the corresponding Dirac equations are almost 
the same which can be seen from formulae (10). Further, it can be shown exactly, that 
Dirac equations in the CJ and MNK versions are equivalent for the states : So and 3 Si 
when mixture of states 3 Si and 3 Di is neglected. 

• For uc and sc bound systems Gross version works better what is related to the large 
difference of constituent masses. 

• The area of the existence of the stable solutions is enlarged with increasing of the con- 
stituent masses which is theoretically understandable. 

• Masses of the bound qq systems obtained from solutions of Dirac equations are bigger then 
masses corresponding to 3D-equations obtained from BS equation for all versions except 
GR version case. 

Note, that for x = 0.5 the stable solutions always exist. 
Acknowledgements. The authors thank A.Rusetsky for useful discussions and comments. 
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Table 1: The dependence of the qq system mass for light constituent quarks on the mixing 
parameter x in the different 3D-reductions of Bethe-Salpeter equations and corresponding Dirac 
equations. "*" denotes the absence of the stable solutions. Masses are given in MeV. 
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Table 2: The dependence of the qq system mass for heavy constituent quarks on the mixing 
parameter x in the different 3D-reductions of Bethe-Salpeter equations and corresponding Dirac 
equations. "*" denotes the absence of the stable solutions. Masses are given in MeV. 
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